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ABSTRACT 

We present the results of a systematic Giant Metrewave Radio Telescope (GMRT) 
survey of 21-cm absorption in a representative and unbiased sample of 35 strong Mg ii 
systems in the redshift range: Zabs~l-10— 1.45, 33 of which have Wr >1 A. The survey 
using '--^400hrs of telescope time has resulted in 9 new 21-cm detections and stringent 
21-cm optical depth upper limits (median 3a optical depth per 10kms~^ of 0.017) 
for the remaining 26 systems. This is by far the largest number of 21-cm detections 
from any single survey of intervening absorbers. Prior to our survey no intervening 
21-cm system was known in the above redshift range and only one system was known 
in the redshift range 0.7 < z < 1.5. We discuss the relation between the detectability 
of 21-cm absorption and various properties of UV absorption lines. We show that if 
Mg II systems are selected with the following criteria, Mg ii doublet ratio <1.3 and 
Wr(Mg i)/Wr(Mg ii) > 0.3, then a detection rate of 21-cm absorption up to 90% can 
be achieved. We estimate n2i, the number per unit redshift of 21-cm absorbers with 
Wr {Mg ii) > Wo and integrated optical depth T21 > %, and show that 7121 decreases 
with increasing redshift. In particular, for Wo = 1.0 A and % > 0.3 km s~^, 7121 falls 
by a factor 4 from < z > — 0.5 to < z > = 1.3. The evolution seems to be stronger 
for stronger Mg 11 systems. Using a subsample of systems for which high frequency 
VLBA images are available, we show that the effect is not related to the structure of 
the background radio sources and is most probably due to the evolution of the cold 
neutral medium filling factor in Mg 11 systems. We find no correlation between the 
velocity spread of the 21-cni absorption feature and Wr(Mg 11) at z 1.3. 
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1 INTRODUCTION 

Observations of high-z galaxies suggest that the global co- 
moving star-formation rate density peaks at 1 < z < 2 and 
then sharply decreases towards 2 ~ (e.g. Madau et al. 
1996, Hopkins 2004). The determination of the mass den- 
sity of the gas and its content (molecules, dust and cold H I 
gas) over the same redshift range provides an independent 
and complementary understanding of the redshift evolution 
of star-formation at similar epochs. While the H I content of 
galaxies can be best probed by surveys of 21-cm emission, 
limited sensitivity of current radio telescopes does not allow 
them to reach beyond the local Universe (Zwaan et al. 2005 
for z ~0; Verheijen et al. 2007 and Catinella et al. 2008 for 
direct detections at z ~0.2; Lah et al. 2007 for statistical 
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detection of H l at 2 ~ 0.24). On the contrary, detection of 
H I in the spectra of distant QSOs in the form of damped 
Lyman- Q absorption provides a luminosity unbiased way of 
probing the evolution of the H I content in the universe 
(Lanzetta et al. 1991; Wolfe et al. 1995; Storrie-Lombardi et 
al. 1996; Ellison et al. 2001; Peroux et al. 2003; Prochaska, 
Herbert-Ford & Wolfe, 2005; Rao, Turnshek & Nester, 2006; 
Prochaska & Wolfe 2008; Noterdaeme et al. 2009). 

A major improvement in the DLA statistics at 2>2 has 
been achieved by the Sloan Digital Sky survey (SDSS). Spec- 
troscopic catalogs are available for DR5 (see Prochaska & 
Wolfe 2008) and for DR7 (see Noterdaeme et al. 2009). For 
lower redshifts, initial detections of DLAs at z < 2 were per- 
formed using the International Ultraviolet Explorer (lUE) 
and the Hubble Space Telescope (HST, Lanzetta, Wolfe & 
Turnshek, 1995; Jannuzi et al. 1998). However, a large num- 
ber of low-z DLAs have been identified after candidate DLAs 
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are pre-selected through the presence of strong metal ab- 
sorption lines (Rao, Turnshek & Nestor, 2006; RTN06 from 
now on). 

It is believed that physical conditions in the neutral 
gas associated with normal galaxies arc influenced by the 
radiative and mechanical feedback from the in-situ star for- 
mation activity (McKee & Ostriker 1977). Thus, it is a good 
idea to investigate whether there is any relationship be- 
tween the evolution of the star-formation rate density and 
the physical conditions in the H I gas probed by DLAs. 
Our understanding of physical conditions in DLAs at z^2 
is largely based on the analysis of H2 and/or atomic fine- 
structure transitions. A systematic search for molecular hy- 
drogen in DLAs at high rcdshifts (zabs^C^LS), using Ultravi- 
olet and Visual Echelle Spectrograph (UVES) at the Very 
Large Telescope (VLT) down to a detection limit of typi- 
cally N{H2) ~ 2 X 10^* cm~^, has resulted in a detection in 
~10-15% of the cases (Ledoux, Petitjcan & Srianand 2003; 
Noterdaeme et al 2008a). HD and CO molecules are detected 
in, respectively, two and one DLAs (Varshalovich et al. 2001; 
Noterdaeme et al. 2008b, Srianand et al. 2008a). The rota- 
tional excitation of the molecules and fine-structure excita- 
tions of C I and C 11 in those systems are consistent with the 
absorption lines originating from dense cold gas in radiative 
equilibrium with the star light produced by the in-situ star- 
formation (Srianand et al. 2000; Ge, Bechtold & Kulkarni 
2001; Reimers et al. 2003; Srianand et al. 2005; Hirashita 
& Ferrara 2005; Cui et al. 2005; Noterdaeme et al. 2007a, 
2007b). Absence of H2 in most of the high-z DLAs is consis- 
tent with the gas being either warm neutral or waxm ionized 
(Petitjean, Srianand & Ledoux 2000; Srianand et al. 2005). 
The absorption lines of C 11* are detected in all the H2 bear- 
ing DLAs (Srianand et al. 2005) and in roughly 50% of the 
DLAs without H2 (Wolfe, Prochaska & Gawiser 2003; Wolfe, 
Gawiser & Prochaska 2003; Srianand et al. 2005; Wolfe et al. 
2008). As C II* is an important coolant, its column density 
combined with (H l) can be used to discuss the cooling rate 
in the absorbing gas and therefore the SFR associated with 
the DLA. Unfortunately for the time being the above men- 
tioned tracers can not be used to probe the physical state 
of the absorbing gas at z^l.S because the useful transitions 
are located below the atmospheric cut-off. 

Most of our initial understanding of the physical state 
of H I gas in the Galactic interstellar medium (ISM) is based 
on the 21-cm absorption line (see Kulkarni & Heiles 1988 for 
details). It is widely believed that the spin temperature (Ts) 
of the H I gas is a reliable tracer of the kinetic temperature 
(see Roy, Chengalur & Srianand 2006 for direct confirma- 
tion). The H I column density for an optically thin cloud 
that covers a fraction /c of the background radio source is 
related to the 21-cm optical depth t{v) in a velocity interval 
v—v+dv and to the spin temperature (Tg) by (e.g. Kulkarni 
& Heiles 1988) 

^■(H I) = 1.835 X 10^* ^ / r{v) dv cm"^ (1) 

/c J 

In addition, the width of the 21-cm line detected in a high 
resolution spectrum yields a direct measurement of (or strin- 
gent upper limit on) the kinetic temperature. Thus detecting 
21-cm absorption is important to probe the physical condi- 
tions in the interstellar medium of galaxies at intermediate 
redshifts. 



RTN06 have shown that DLAs essentially have Mg 11 
rest equivalent width, ^^(Mg llA2796)>0.6A. Therefore, the 
search of 21-cm absorption in a sample of strong Mg II ab- 
sorbers is an unique way to probe the redshift evolution of 
physical conditions in DLAs like absorption systems at low- 
z. In the next Section we summarize the previous efforts to 
search for 21-cm absorption in DLAs and Mg 11 systems. 
Due to vaxious practical reasons (e.g. frequency coverage 
of receivers, radio frequency interference, RFI) the redshift 
coverage is sparse and measurements are available only for a 
few systems, especially in the redshift range: 1 < 2: < 2. The 
Giant Meterwave Radio Telescope (GMRT), with its very 
sensitive 610 MHz receiver, is the only radio telescope avail- 
able at present in the relatively RFI-clean environment (e.g. 
compared to Green Bank Telescope or Westerbork Synthesis 
Radio Telescope) for covering part of this redshift range. In 
addition, the large number of publicly available high- 2; QSO 
spectra from SDSS allow one to construct the large samples 
of Mg II absorbers suitable for the searches of DLAs/21-cm 
absorbers (cf. Section 3). Motivated by this we have con- 
ducted a systematic search for 21-cm absorption in a rep- 
resentative sample of Mg 11 systems in the redshift range 
1.10 < z < 1.45 and discovered 9 new 21-cm absorbers. 
This is by far the highest number of detections from any 
single survey of intervening 21-cm absorption (see Srianand, 
Gupta & Petitjean 2007; Gupta et al. 2007 and Srianand et 
al. 2008b for preliminary results from the survey). 

This paper is structured as follows. We give a brief sum- 
mary of the past 21-cm searches performed by other groups 
in Section 2. In Section 3, we describe our sample and details 
of observations are presented in Section 4. In Section 5, we 
describe in detail the individual 21-cm absorbers detected 
in our survey. Detectability of 21-cm absorbers and its rela- 
tionship to metal absorption equivalent widths, equivalent 
width ratios and structure of QSOs at radio wavelengths are 
discussed in Section 6. In Section 7, we derive the number 
density per unit redshift of 21-cm absorbers. Using the low-2 
sample of Lane (2000) we investigate the redshift evolution 
of the number density of 21-cm absorbers. In Section 8, we 
explore the relationship between the velocity spread of 21- 
cm absorption and the equivalent width of the Mg 11 absorp- 
tion line. Discussion of our results along with a summary are 
presented in Section 9. 



2 BRIEF SUMMARY OF PREVIOUS 21-CM 
ABSORPTION LINE SURVEYS 

There have been several systematic surveys for 21-cm ab- 
sorption in DLAs and/or strong metal line absorption 
systems undertaken by various groups over the past two 
decades. Briggs & Wolfe (1983) searched 18 Mg 11 systems, 
with W^r(Mg ii)>0.5A and 0.3 < z < 1.8, and detected 
two. Lane (2000), in her thesis work, searched 62 Mg 11 
systems at 0.2 < z < 0.9 without imposing any equiva- 
lent widths cut-off and detected 21-cm absorption in 3 sys- 
tems. Results of 21-cm searches are available for 7 DLAs at 
z > 2.8 (Carilli et al. 1996; Kanekar & Chengalur 2003) 
with only one detection (Kanekar et al. 2007). Measure- 
ments are consistent with Ts > 1000 K. There are five 21- 
cm detections reported in the literature for 1.7 < z < 2.8 
[zabB= 2.289 towards TXS0311-h430 (York et al. 2007); 
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Figure 1. Redshift distribution of Mg II systems that were 
searched for 21-cm absorption. The filled histogram is the GMRT 
sample of 35 Mg II systems presented in this paper (33 of these 
absorption systems have Wi{Mg llA2796)^lA). The solid line his- 
togram is for the sample of Lane (2000). The hatched histogram 
corresponds to 21-cm detections in this sample. The distribution 
for the Wi{Mg llA2796)> lA sub-set of these systems is given by 
the dashed line histogram. 
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Figure 2. Top panels: The Mg II doublet ratios and Mg IIA2796 
equivalent widths in our sample are plotted against the 1.4 GHz 
peak flux density of the background source. Triangles are for sys- 
tems from DR3 (P06) and circles are for additional systems in 
DR5. Filled and open symbols are for systems with and with- 
out GMRT observations. Typical error in the measurements are 
shown in these panels. Third and bottom panel: Distribution of 
radio flux densities in our sample. The shaded histograms gives 
the distributions of objects observed with GMRT. 



«abs= 2.347 towards B0438-436 (Kanekar et al. 2006); Zabs = 
2.039 towards B0458-020 (Wolfe et al. 1985); Zabs= 1-943 to- 
wards 1157+014 (Wolfe & Briggs, 1981); and Zabs = 1.776 to- 
wards 1331-1-170 (Wolfe & Davis, 1979)]. But non-detections 
are not systematically reported in this redshift range. At 
z < 0.9, Kanekar & Chengalur (2003) summarize 21-cm op- 
tical depths for 14 known DLAs and 12 detections (including 
their own observations of two DLAs namely PKS 1629-f 12 
and PKS 2128-13). Curran et al. (2007) report the detection 
of a 21-cm absorption at 2abs= 0.656 in the complex sight- 
line towards 3C336. In total 18 intervening 21-cm absorption 
systems are known till now. 

We plot in Fig. [T]the redshift distribution of the Mg ll 
absorption systems surveyed by Lane (2000). This is the 
only large survey at low-z for which both detections and 
non-detections are systematically reported. It includes 62 
observations using WSRT and 10 other systems from the 
literature satisfying their selection criterion (see Lane 2000 
for details). The detections shown as a hatched histogram 
include detections reported in Lane (2000) together with de- 
tections from better quality data by Kanekar & Chengalur 
(2003) and Curran et al. (2007) for systems that were orig- 
inally reported as non-detections. In the same figure, the 
filled histogram shows the distribution of Mg ll systems in 
our GMRT sample. For equivalent width cutoff of ~ lA, our 
GMRT sample has more than twice the number of systems 
investigated by Lane (2000). A detailed comparison of the 
results from low and high redshift samples will be discussed 
in Section 7. 



3 SAMPLE OF Mg ii SYSTEMS IN FRONT OF 
RADIO QUASARS 

RTN06 have shown that DLAs can be preselected on the 
basis of the equivalent widths of Mg ll, Fe ll and Mg I ab- 
sorption lines. Specifically, they found that 36%±6% of the 
Mg II absorbers with Fe iiA2600 and Mg nA2796 rest equiv- 
alent widths greater than 0.5 A are DLAs. Motivated by 
this we have conducted a GMRT 21-cm absorption survey 
over the past two years in absorption systems selected on 
the basis of Mg llA2796 rest equivalent width in the redshift 
range 1.10< z < 1.45. The constraint on the redshift range 
comes from the frequency range covered by the GMRT 610 
MHz receiver. Our sample is drawn from the identification 
of Mg 11 systems by Prochter, Prochaska & Buries (2006, 
hereafter P06) in SDSS DR3 and by us using our automatic 
procedure for additional systems in DR5. We do not make 
any further selection cut based on the equivalent widths of 
Fe II and Mg I absorption lines. Thus our candidate selection 
is based on Wr{Mg llA2796) alone. 

There are 1953 Mg ll systems in P06 with 
Wr{Mg IIA2796) > 1.0 A in the redshift range: 1.10< z < 
1.45. By cross-correlating the corresponding QSO positions 
with the NRAO VLA Sky Survey (NVSS) radio catalog 
(Condon et al. 1998) we found 45 of these Mg ii systems 
having background QSO with flux density at 1.4 GHz in ex- 
cess of 50 mjy. We manually checked the SDSS spectra of all 
these sources to confirm the Mg ll systems and remeasure 
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accurately the A2796 rest equivalent width as well as other 
lines. This step is essential as the identifications of Mg ll sys- 
tems are not 100% complete and errors in the Mg ll equiva- 
lent width measurement can be large (see discussion on this 
in P06) when one is using automatic procedures. We find 
that Mg II identification in one of the QSOs (i.e 2abs= 1.109 
towards the very bright radio source J104833. 69-1-600845. 7) 
is spurious. While manually checking the SDSS spectra we 
detected an additional (i.e. not catalogued in POO) Mg ll sys- 
tem at Zabs = 1.3739 towards J123431. 72-1-645556.5. Since 
this system has Wr(Mg iiA2796) = 0.95±0.07, we included 
it in our sample as it satisfies the selection criteria of 
Wr>lk within 1 a and 1.10<2abs<1.45. 

In order to minimize the possibility that the opti- 
cal and radio sightlines are not identical we looked at 
the ~5" resolution radio images of the 44 remaining 
sources (with 45 absorption line systems) obtained in the 
Faint Images of the Radio Sky at Twenty-centimeters 
(FIRST) survejQ (White et al. 1997). Eight sources in 
the list (J024534.05-h010813.8, J095631.05-f-404628.2, 
J095908.31~f 024309.6, J103838.82+494736.8, 
J130907.98+522437.3, J141701.94-H541340.6, 
J160846.76-h374850.6 and J223334.89-075043.2) are re- 
solved in the FIRST survey maps and the associated radio 
component closest to the optical source has a flux density 
less than 50 mjy. For this reason these sources are not in- 
cluded in our sample. In the case of J124538. 35-1-551132. 66 
and J144553. 46-1-034732. 4 the radio peak of interest is more 
than 8" away from the corresponding optical positions. 
As these separations are much larger than the astrometric 
accuracy of FIRST and SDSS we do not include these 
sources in our list. Another source, J100121. 07-1-555355. 8, 
is a gravitational lens commonly known as Q0957-|-561A,B. 
The two optical images are separated by ~ 6". The radio 
structure is complex with core and jets. At 1420 MHz, 
the radio image is deconvolved into three components and 
none of the peaks coincide with any of the two optical 
images. The dominant radio component is located more 
than 4" from the optical sources. For these reasons we did 
not include this source in our sample. Note that 21-cm 
absorption associated with the Zabs = 1.3911 Mg ll system 
has been searched for at GMRT by Kanekar & Chengalur 
(2003). 

We ended up with a complete sample of 34 Mg ll 
(Wr>lA) systems (in front of 33 radio sources) with flux 
density above 50mJy from SDSS DR3. In our survey we 
have searched for 21-cm absorption in 24 of these systems. 
The distribution of the radio flux density in these sources 
is given in the third panel of Fig. [5] From this figure it is 
clear that we have searched for 21-cm absorption in all but 
the two systems with background sources having fiux den- 
sity greater than 100 mJy. We are only ~45% complete for 
fiux densities in the range 50—100 mJy. 

Using an automatic procedure (similar to the one used 
by P06), we have identified 1548 additional strong Mg ll sys- 
tems in the same redshift range from DR5 of SDSS. Follow- 
ing the steps discussed above we have identified 29 systems 



^ FIRST images are not available for the two sources, namely, 
J021452.29-I-140527.4 and J123431. 72-1-645556.5. We retain these 
in our sample. 



in front of 29 QSOs having at least a Si.4ghz > 50 mJy com- 
ponent coinciding with the optical position in the FIRST 
image. Due to scheduling constraints, we could only cover 8 
of these sources in our GMRT observations. 

Therefore in total we have observed a sample of 32 
Wv^lK Mg II systems in front of the 31 radio sources drawn 
from SDSS. In the bottom panel of Fig [2] we give the dis- 
tribution of fiuxes in our observed sample. In summary, we 
have observed, respectively, ~70% (resp. 35%) of the avail- 
able sources with fiux densities 5'i.4ghz > 100 mJy (resp. in 
the range 50-100 mJy). From the top two panels of Fig[2]it is 
apparent that the distribution of Wr(Mg llA2796) and Mg ll 
doublet ratio (DR) in the systems observed with GMRT 
(filled circles) represent the parent distribution well. The 
two distribution Kolmogorov-Smirnov test (KS test) gives 
the probability of 70% and 36% respectively for the dis- 
tribution of M^r(Mg IIA2796) and DR(Mg ll) in the sources 
with GMRT observations to represent the distribution of the 
parent population. Therefore, despite incompleteness in our 
survey towards low flux density levels (i.e <100 mJy), the 
sample is representative of the parent population as far as 
Wi{Mg IIA2796) and DR(Mg ll) are concerned. 

In addition, we have observed three more sources i.e. 
J010826.84-003724.1 (York et al. 2006), and J0240-2309 
and J160456. 14-001907.1 (Lanzetta et al. 1987) because 
they have a total fiux density at 610 MHz in excess of 1 
Jy and are associated with strong Mg ll absorptions. We 
notice that addition of these systems (that were observed in 
the early stages of our survey) does not change the statis- 
tical conclusions drawn above. Thus, in total we have ob- 
served 35 Mg II systems, 33 of which have Wr(Mg ll)^ lA, 
to search for 21-cm absorption. In Table[l]we give the list of 
observed sources together with the emission redshift, Zcm in 
column 2, the absorption redshift (zabs) in column 3, the 
rest equivalent widths of Mg llAA2796,2803 (columns 4 and 
5), Mg IA2852 (column 6) and Fe llA2600 (column 7), the ra- 
dio fiux density at 1420 MHz (column 8), morphology of the 
radio sources (column 9) and the origin of the Mg ll systems 
(column 10). The equivalent width estimates in Table [T] were 
obtained by integrating the absorption over the absorption 
profile. For consistency these measurements were compared 
against estimates from Voigt profile fits. Differences of the 
order of 0.2A or less were found between estimates from 
these methods. This is similar to uncertainties (~0.3A) from 
estimating equivalent widths by automatic procedures (see 
P06). 

Preliminary results from the first phase of the GMRT 
survey reporting 3 new detections in 10 systems were pre- 
sented in Gupta et al. 2007 (see also Srianand, Gupta & 
Petitjean, 2007 for details regarding J0240-2309). Detailed 
analysis of two detections from our survey with strong sig- 
natures of dust reddening is presented in Srianand et al. 
(2008b). Here we present the results for the entire survey. 



4 GMRT OBSERVATIONS 

We observed 35 Mg ll systems with GMRT 610-MHz band 
using in total ~400 hrs of telescope time. Details of the ob- 
servations are given in Tabled For our survey, we have usu- 
ally used the 1 MHz baseband bandwidth split into 128 fre- 
quency channels yielding a spectral resolution of ~4kms~^. 



Intermediate redshift 21-cm absorbers 5 



Table 1. Sample of Mg II systems observed with GMRT. 



Source name 


•2cm 




M/r(Mg IIA2796) 


VVr(Mg IIA2803) 


VVr(Mg IA2852) 


Wi(Fc nA2600) 


Flux 


Morp. 


Ref. 










(A) 


(A) 


/ X ^ 

(A) 


(A) 


(mJy) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


J010826. 84-003724.1 


1.373 


1 


.3710 


0.43±0.05 


0.39±0.05 


0.20±0.05 


0.33±0.06 


894.9 


c 


Y 


J015454. 36-000723. 2 


1.828 


1 


.1803 


1.37±0.19 


1.23±0.19 


<0.22 


0.65±0.26 


245.8 


C 


P 


J021452.29-I-140527.4 


2.178 


1 


.4463 


2.53±0.22 


1.55±0.19 


1.04±0.23 


1.28±0.23 


119.5 


CN 


P 


J0240-2309t 


2.223 


1 


.3647 


1.85±0.002 


1.65±0.002 


0.27±0.002 


0.87±0.003 


6027 


C* 


L 


J025928.51-001959.9 


2.000 


1 


.3370 


1.76±0.04 


1.35±0.05 


<0.05 


0.36±0.05 


226.7 


c 


P 


J074237.38-I-394435.6 


2.200 


1 


.1485 


2.48±0.41 


1.41±0.35 


<0.26 


<0.69 


99.3 


c 


P 


J074809. 46-1-300630. 5 


1.695 


1 


.4470 


3.68±0.06 


3.28±0.07 


0.85±0.08 


2.17±0.08 


200.4 


c 


P 


J080248.43-I-291734.1 


2.380 


1 


.3648 


1.04±0.19 


0.70±0.16 


<0.79 


<0.23 


190.4 


c 


P 


J080442. 23-1-301237.0 


1.452 


1 


.1908 


1.28±0.05 


1.07±0.06 


0.19±0.06 


1.05±0.07 


1040 


R 


P 


J080839. 66-1-495036. 5 


1.436 


1 


.4071 


1.33±0.09 


1.24±0.08 


0.48±0.12 


0.60±0.12 


846.8 


C 


T 


J081247.78-I-252242.0 


1.803 


1 


.3683 


3.18±0.10 


3.11±0.11 


0.87±0.12 


2.22±0.14 


69.6 


C? 


P 


J084506. 24-1-425718. 3 


2.095 


1 


.1147 


1.24±0.09 


1.34±0.09 


<0.11 


0.87±0.09 


203.4 


C 


P 


J085042.24-I-515911.7 


1.894 


1 


.3265 


4.56±0.12 


4.02±0.12 


2.19±0.12 


3.44±0.15 


61.2 


C 


P 


J085244.74-I-343540.4 


1.655 


1 


.3095 


2.91±0.11 


2.83±0.11 


1.27±0.12 


2.26±0.18 


66.7 


C 


T 


J095327.95-I-322551.6 


1.575 


1 


.2372 


1.55±0.05 


1.49±0.05 


0.35±0.06 


1.18±0.08 


127.4 


C 


T 


J100842.71-I-621955.7 


1.875 


1 


.2080 


4.04±0.22 


3.03±0.20 


1.37±0.30 


2.38±0.28 


73.8 


R 


P 


J101742. 63-1-535635.0* 


1.397 


1 


.3055 


2.68±0.10 


2.62±0.11 


0.70±0.11 


1.98±0.14 


64.0 


R 


P 


J102258.41-I-123429.7 


1.727 


1 


.2505 


3.06±0.07 


2.77±0.08 
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Col. 1: SDSS source name; col. 2: QSO emission redshift; col. 3: Absorption redshift of Mg II systems as determined from 
Mg llAA2796,2803, Mg lA2852and Fe nA2600absorption lines; cols. 4, 5, 6, 7: Rest equivalent widths or la limit of Mg IIA2796, 
Mg IIA2803, Mg IA2852, Fe llA2600 absorption lines respectively; col. 8: Peak flux density at 1.4 GHz as measured from FIRST (or 
NVSS when FIRST image is not available) survey, col. 9: Radio morphology as determined from FIRST (or NVSS) where C=compact 
at the FIRST resolution (~5"), CN=compact at NVSS resolution (~45"), R=resolved in the FIRST image, D=double lobed and 
T=triple radio source, and col. 10: Origin of the Mg II systems, where Y is for York et al. (2006), P is for P06, L for Lanzetta et al. 
(1987), T for 'this paper' i.e. drawn from SDSS using a procedure developed by one of us and M is for detected while manually 
checking the SDSS spectra for other systems. 

^ Absorption line equivalent widths for this system are estimated from the VLT UVES data presented in Srianand et al. (2007). 
t Source is not covered by the FIRST and NVSS surveys but is a VLA phase calibrator for A, B, C and D conflgurations. 
* This source is extended in FIRST but compact in our GMRT image. 



Broader (2 MHz for J0804-I-3012) and narrower (0.5 MHz 
for J2358-1020 and 0.25 MHz for J2340-0053) bandwidths 
were used as well for confirming the detected 21-cm absorp- 
tions. During each observations, the local oscillator chain 
was tuned to centre the baseband at the redshifted 21-cm 
frequency defined by the mean redshift estimated from the 
Mg 11 and other absorption lines detected in the SDSS spec- 
trum. Data were acquired in the two circular polarization 
channels called RR and LL. We typically observed each 
absorption line system for 8hrs which included calibration 



overheads. For flux and bandpass calibration we observed 
one of the standard flux density calibrators 3C 48, 3C 147 
and 3C286 for 10-15 mins every two hours. A phase cali- 
brator was also observed for 10 mins every ~40 mins to get 
reliable phase solutions. Total on-source time after exclud- 
ing the telescope set-up time and calibration overheads are 
provided in Table (2] 

GMRT data were reduced using the NRAO Astronomi- 
cal Image Processing System (AIPS) following the standard 
procedures as described in Gupta et al. (2006). After the 
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Figure 3. GMRT maps of sources resolved in our observations with contour levels as nx (-1,1,2,4, ...) mjy beam ^. J0804+3012: 0.84 
mjy beam~^ rms and beam of 6.36" x4. 09", position angle PA=67° and n=4.0. J0812+2522: 0.90 mJy beam'^ rms and beam of 
5.60"x4.37", PA=66° and n=5.0. J1008+5219: 0.51 mjy beam"! rms and beam of 8.87"x4.92", P.A.=67° and n=1.6. J1022+1234: 
0.90 mJy beam"! rms and beam of 6.13"x4.49", PA=45° and n=6.0. J1145+0455: 0.90 mjy beam-^ rms and beam of 5.85"x4.70", 
PA=63° and n=3.0. J1329+1053: 0.70 mjy beam"! rms and beam of 6.85" x5. 88", PA=-26° and n=3.2. J1456+0005: 0.48 mjy beam'^ 
rms and beam of 5.93"x4.78", PA=59° and n=2.0. J1510+5958: 1.10 mjy beam"! rms and beam of 5.67"x4.71", PA=36° and n=6.0. 
J2312-0109: 0.95 mjy beam"! rms and beam of 7.38"x4.66", PA=53° and n=4.0. 



initial flagging and calibration, source and calibrator data 
were examined to flag and exclude the baselines and times- 
tamps afltected by Radio Frequency Interference (RFI). Ap- 
plying complex gains and bandpass obtained using the flux 
and phase calibrators, a continuum map of the source was 
made using absorption-free channels. Using this map as a 
model, self-calibration complex gains were determined and 
applied to all the frequency channels. The same continuum 
map was then used to subtract the continuum emission from 
the dataset. This continuum subtracted dataset was then 
imaged separately in stokes RR and LL to get 3-dimensional 
(with third axis as frequency) data cubes. The spectra at the 
quasar positions were extracted from these cubes and com- 



pared for consistency. If necessary a first-order cubic-spline 
was fitted to remove the residual continuum from spectra. 
The two polarization channels were then combined to get 
the final stokes I spectrum, which was then shifted to helio- 
centric frame. 

We have obtained useful radio spectra for 35 Mg ll sys- 
tems in front of 34 radio sources. Of these 34 radio sources 
20 are compact in FIRST (deconvolved sizes < 2") and 
our 610 MHz images. These are assigned morphology C 
in Tables [1] and [3] Of the remaining 14 sources, for two 
sources (i.e. J0214-hl405 and J1234-H6455) FIRST images 
are not available. Both sources are compact in NVSS and 
our GMRT 610-MHz images. Radio emission associated with 
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Table 3. Observational results for Mg II systems observed with GMRT. 
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.754 
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Col. 1: Source name; col. 2: Mg II absorption redshift; col. 3: Morphology of the radio source as determined from our observations; col. 
4: Peak flux density in mJy; col. 5 and 6: Spectral resolution in kms~^ and spectral rms in mJy beam~^ channel"^; col. 7: Maximum 
of optical depth in one channel or Icr limit to it; col. 8: 3cr optical depth limit for spectra smoothed to ~10kms~^; col. 9: integrated 
21-cm optical depth or 3ct upper limit in case of non-detections considering (5?;=10 km s~^. 



J0214-fl405 in our GMRT image (with beam 5.5" x4.6" and 
PA=57°) is well represented by a Gaussian component im- 
plying a deconvolved size of 2.5" x 1.6" with position angle 
56° . Despite the deconvolved radio source size being slightly 
greater than 2", >90% of the total emission is contained 
in the above mentioned Gaussian component, so we con- 
sider this source to be compact. Gaussian deconvolution of 
J1234-I-6455 yields a size of 0.7"x0.0" with position angle 
12° implying that the source is compact. This source is also 
compact in sub-arcsecond scale CLASS images at 8.4 GHz. 
Thus there are 22 sources compact at arcsec scales in our 
sample and 23 Mg ll systems along their lines of sight with 
GMRT data. From Table |3] one can see eight (i.e ~35%) of 
these systems show detectable 21-cm absorption. 

Eleven sources in our sample show extensions in our 
GMRT continuum maps (see Fig. [3] in this article and Fig. 
4 in Gupta et al. 2007). For 8 of these sources spectra were 



extracted at the position of the peak flux that also coin- 
cides well with the optical position. For the three remaining 
sources (J1604-0019, J1411-0300 and J1510+5958) we also 
extracted spectra from the hot spots that are resolved from 
the core in our GMRT observations. In none of these cases 
we have detected 21-cm absorption towards hot spots. We 
address the issue of radio morphology of our target sources 
and its influence on the detectability of 21-cm absorption in 
following sections. 

Absorption profiles of the nine newly discovered 21-cm 
absorbers are presented in Fig. |4] whereas spectra with no 
detectable 21-cm absorption are shown in Fig.[S] Results are 
summarised in Table |3] In the 3rd to 6th columns of this ta- 
ble we provide the morphology, the peak flux {Fp) measured 
in our GMRT 610 MHz continuum image, the velocity res- 
olution (Sv) and the measured spectral rms. In the seventh 
column, we provide the maximum optical depth per channel 
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Figure 4. GMRT spectra of detected 21-cm absorption lines. Individual Gaussian components and resultant fits to the absorption 
profiles are overplotted as dotted and continuous lines respectively. Residuals, on an offset arbitrarily sfiifted for clarity, are also shown. 



or the 1(7 limit on it. The next column gives the 3(t optical 
depth limit obtained as follows: rso- = —ln{l — 3AF/Fp), 
where Fp is the peak flux and AF is the spectral rms for the 
spectrum smoothed to 10 km s~^resolution. This smoothing 
velocity scale is chosen to enable comparative studies with 
the existing low redshift (i.e Zabs<l-0) measurements (see 
Section 7). In the case of detections, we measure this quan- 
tity in the line free spectral range. The last column in the 
table gives the integrated optical depth or the limit on it. In 
the later case we used a Gaussian width of 10 km s~^ for 
the integration (i.e 10.6 xrso-). 

For detections, it is important to rule out the possibil- 
ity of false identification even in case of good redshift coin- 
cidence. For spectra with high SNR (i.e 21-cm absorption 
lines towards J0108-0037, J0850-F5159 and J1058+4939) 
we split the data into different time and UV ranges and 



check for consistency of the spectra. We also compared the 
spectra with spectra of other bright sources in the field. For 
spectra of poor SNR or when the narrow absorption line fea- 
ture is spread over only one or two frequency channels (i.e. 
along the lines of sight towards J0804+3012, J0808+4950, 
J0852+3435, J1623+0718, J2340-0053 and J2358-1020), 
observations were repeated at a different epoch. In all cases 
we observed a consistent absorption profile at a frequency 
that is in good agreement with the Doppler shift expected 
from the heliocentric motion. For J0808+4950, J0852+3435 
and J1623-I-0718 spectra at different epochs were obtained at 
the same spectral resolution 4 km s~^) and the final spec- 
tra presented in Fig. |3] are weighted averages of these. The 
absorption line system towards J0804+3012 was observed 
three times. The later two observations were done using 
2 MHz bandwidth split into 128 channels to obtain a wider 
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Figure 5. GMRT spectra of sources with non-detection. Arrows mark the expected positions of 21-cm absorption hnes based on metal 
absorption lines. The shaded regions mark frequency ranges affected by RFI. The 2:abs=l-3739 and 1.3829 absorbers towards J1234+6455 
are referred to as A and B respectively. For J1411-0300, PI and P2 correspond to spectrum towards two peaks (see Gupta et al. 2007 
for detail). The spectra toward northern and southern radio peaks of J1510+5958 are marked as N and S. 



spectral coverage. The final spectrum presented here was ob- 
tained by smoothing the spectrum from the first observing 
run and combining this with the later two spectra of spec- 
tral resolution ~7km s~^(see also Gupta et al. 2007). Nar- 
row 21-cm absorption profiles detected towards J2340— 0053 
and J2358— 1020 were reobserved at higher spectral resolu- 
tion. Here, we present only the highest resolution spectra. 



5 DETAILS OF NEW 21-CM ABSORBERS 

In this Section, we give details on the nine 21-cm absorbers 
detected in our survey (Fig. |4]). The total integrated 21-cm 
optical depths and velocity spreads of the absorption (AV^) 



are listed in Table 21 We compute AV using the apparent 
optical depth profile (see Fig. |6]l and the method described 
by Ledoux et al. (2006) . It is usual procedure to decompose 
the 21-cm absorption profile into multiple Gaussian func- 
tions. For a Gaussian profile A''(H l) and the peak optical 
depth (rp) are related by. 



N(m) = 1.93 X 10^® Tp ^ At; cm"^ 



(2) 



Here, Tp , Av, and /c are the peak optical depth, the FWHM 
in km s~^of the fitted Gaussian, and the covering factor of 
the absorbing cloud respectively. Individual Gaussian fits 
and fit parameters are presented in Fig. [4] and Table [4] re- 
spectively. 
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Figure 5. Continued. 



In principle one can constrain the kinetic temperature 
Tk using the width of the absorption lines (Tk < 21.86 x 
Av^K; e.g. Heiles & Troland 2003). However, FWHMs of 
individual Gaussian components fitted to the absorption 
profiles are generally broad and range from 5 km s~^to 
100 km s~^ when our resolution is ~4 km s~^. Most proba- 
bly these line widths are dominated by processes other than 
simple thermal broadening and constraints on Tk are poor. 
J0108-0037: The 21-cm absorption profile is weU fitted 
with a single Gaussian function with rp=0.070±0.002 and 
An = 16±1 km s~^. The absorption line width corresponds 
to an upper limit on the kinetic temperature of Tk < 5600 K. 
Radio source J0108— 0037 is resolved into two components 
separated by ~5mas (~40pc at Zabs = 1.3710) in the VLBA 
images at 2 and 8 GHz (Beasley et al. 2002). Of these two 
the strongest component contributes about 70% of the total 
fiux. Since it is difficult to ascertain, which of these compo- 



nents is the radio core, estimate of /c is highly uncertain. 
From the integrated optical depth given in Tabled we de- 
rive A^(H l)=2.4xl0^° cm"^ for /c = 1 and = 100 K (as 
in cold neutral medium (CNM) of the Galaxy (Kulkarni & 
Heiles 1988) ). This is a lower limit on N{i{ l) as /c is most 
probably less than 1 and Ts could be larger than 100 K. Thus 
A'^(H l) in the 21-cm component is consistent with this sys- 
tem being a DLA. The Zn ii+Crii blend at A = 2062 A, the 
Si II and the Mn ll absorption lines are clearly detected even 
in the low dispersion spectrum (see Fig. [Tjl . As can be seen 
from Table [H this system has Wr(Mg llA2796)< 0.5 A and 
was observed in the early stages of the survey mainly be- 
cause of the presence of strong absorption lines from weak 
metal transitions (see Fig. [7]). For this reason this system 
will not figure in any analysis we perform using equivalent 
width limited samples. 
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Figure 5. Continued. 



Table 4. Characteristics of 21-cm absorption profiles and multiple Gaussian fits. 



QSO 


J' rdv 


AV 


/cJV(HI)" 


^abs 


Ad'' 




/ciVCHI)" 


J0108-0037 


1.29±0.04 


27.31 


2.37±0.07 


1.37099 


16±1 


0.070±0.002 


0.22±0.02 


J0804+3012 


0.39±0.03 


130.1 


0.72±0.06 


1.19069 


86±47 


0.002±0.001 


0.03±0.02 










1.19109 


43±12 


0.004±0.002 


0.03±0.02 


J0808+4950 


0.10±0.02 


23.82 


0.18±0.04 


1.40732 


11±2 


O.OOSitO.OOl 


0.02±0.01 


J0850+5159 


15.3±0.5 


61.38 


28.1±0.9 


1.32674 


49±16 


0.117±0.036 


1.11±0.50 










1.32692 


24±4 


0.347±0.108 


1.61±0.57 


J0852+3435 


6.91±0.34 


76.13 


12.7±0.6 


1.30919 


23±7 


0.078±0.025 


0.35±0.15 










1.30945 


63±12 


0.079±0.009 


0.96±0.21 


J1058+4939 


0.41±0.04 


32.85 


0.75±0.07 


1.21168 


26±9 


0.006±0.001 


0.03±0.01 










1.21181 


15±2 


O.OlSitO.OOl 


0.05±0.01 


J1623+0718 


o.giito.io 


34.57 


1.67±0.18 


1.33567 


11±2 


0.042±0.005 


0.09±0.02 










1.33585 


17±5 


0.023±0.005 


0.08±0.03 


J2340-0053 


2.50±0.24 


7.76 


4.59±0.44 


1.36087 


2.5±0.5 


0.662±0.207 


0.32±0.12 










1.36090 


5.8±2.8 


0.183±0.042 


0.21±0.11 


J2358-1020 


0.23±0.04 


12.48 


0.42±0.07 


1.17304 


6±1 


0.034±0.004 


0.04±0.01 


in units of 10^** 


cm-^; " FWHM 


in km s~^ ; 


peak optical depth; 


in units of 10 


" cm 
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Figure 7. Normalised rest frame SDSS spectra of the Mg ll systems with 21-cm detection. Vertical dashed lines mark the location of 
interesting metal line transitions. 



J0804+3012: Broad and shallow 21-cm absorption de- 
tected towards this QSO was confirmed by observations at 
different epochs and using a wider baseband bandwidth of 
2 MHz. Therefore, latter spectra have a factor 2 coarser 
velocity resolution (see Table [2)|. In Fig. |4l we show the 
weighted average of the spectra obtained after smoothing the 
spectrum of higher resolution. This system has the broadest 
velocity profile among the 21-cm absorbers in our sample. 
The overall absorption profile is fitted with a main com- 
ponent with Av ~43 km s~^and a broad wing, At; ~86 
km s~^. Based on the integrated optical depth we derive 
iV(H i)=7.0xl0i''(Ts//c)cm-2. 

This source is resolved into three components in the 
MERLIN images obtained at 1.6 GHz and 5 GHz with an- 
gular resolutions 0.178" xO. 155" and 0.073" xO. 037" respec- 
tively (Kunert et al. 2002) . The angular extent of the radio 



emission in these images is ~ 1.4" with 77% of the fiux con- 
tributed by the central component. However, at 1.4 GHz 
these three components account for only 50% of the flux 
measured in FIRST images. The structure of the source is 
thus highly complex. This could very well explain the large 
velocity spread of the 21-cm absorption. This source is ideal 
for VLBA spectroscopy to study the internal structure of 
the intervening Mg ll system. Apart from Mg ll, Fe ll and 
Mg I absorption lines no other metal absorption lines are 
detected in the SDSS spectrum (see Fig. [7|. 

J0808-I-4950: The background QSO at Zcm = 1.432 is clas- 
sified as highly polarized Blazar/QSO with radio compo- 
nents at mas scales showing superluminal motion (see Finer 
et al. 2007). The measured apparent velocity is ~ 15c 
with c being the speed of light. This source is also known to 
show strong fiux variability at high radio frequencies (Teras- 
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Table 2. Observing log for the GMRT observations. 



J010e26.a4-003724.1 J080442.23+301337.0 
0.005 



J0B0a39.6e 4-495036.5 



Source name 




Date 




Time 


On. Width 










(hr) 


/i — 1\ 
(kms ) 


(1) 


(2) 


(3) 




(4) 


(5) 


J0108-0037 


1.3710 


2007 


Jul 03 


4.4 


3.9 


J0154-0007 


1.1803 


2007 


Jan 30 


5.3 


3.6 


J0214+1405 


1.4463 


2006 


Jul 16 


4.4 


4.0 


J0240-2309 


1.3647 


2004 


Feb 09 


8.0 


3.9 


J0259-0019 


1.3370 


2006 


Oct 06 


2.3 


3.9 


J0742+3944 


1.1485 


2007 


Dec 31 


6.3 


3.5 


J0748+3006 


1.4470 


2006 


Aug 01 


5.5 


4.0 


J0802+2917 


1.3648 


2007 


Nov 26 


7.0 


3.9 


J0804+3012 


1.1908 


2006 


Aug 22 


6.0 


3.6 






2006 


Sep 15,16 


6.2 


y 22.128 


J0808+4950 


1.4071 


2006 


Dec 06 


5.3 


4.0 






2008 


May 08 


5.1 


4.0 


J0812+2522 


1.3683 


2007 


Feb 01 


5.8 


3.9 


J0845+4257 


1.1147 


2006 


Sep 09 


5.0 


3.5 


J0850+5159 


1.3265 


2007 


Nov 06 


7.3 


3.8 


J0852+3435 


1.3095 


2007 


Nov 05 


7.5 


3.8 






2007 


Nov 30 


6.9 


3.8 






2008 


Mar 08 


6.2 


3.8 


J0953+3225 


1.2372 


2007 


Jun 05 


5.8 


3.7 


J1008+6219 


1.2080 


2006 


Oct 06 


6.5 


3.6 


J1017+5356 


1.3055 


2006 


Aug 29 


4.4 


3.8 


J1022+1234 


1.2505 


2007 


Nov 29 


8.3 


3.7 


J1058+4939 


1.2120 


2008 


May 10 


5.9 


3.6 


J1126+4516 


1.3022 


2007 


Dec 29 


5.5 


3.8 


Jl 145+0455 


1.3433 


2007 


Jan 30 


5.2 


3.9 


J1208+5441 


1.2110 


2008 


May 11 


6.1 


3.7 


J1232+5722 


1.3429 


2007 


Dec 28 


5.6 


3.9 


J1234+6455 


1.3739 


2008 


Jan 02 


5.6 


3.9 




1.3829 


2007 


Dec 30 


5.4 


3.9 


J1329+1053 


1.1645 


2007 


Jun 21 


5.2 


3.6 


J1411-0300 


1.4160 


2006 


Jul 06,07 


3.9 


4.0 


J1456+0005 


1.3512 


2008 


Jan 01 


6.1 


3.9 


J1508+3347 


1.1650 


2007 


Jun 05 


4.1 


3.6 


J1510+5958 


1.3720 


2007 


Feb 02 


5.6 


3.9 


J1604-0019 


1.3245 


2006 


Jun 19 


3.3 


3g2,256 


J1623+0718 


1.3350 


2007 


Dec 16 


5.4 


3.9 






2008 


May 09 


6.3 


3.9 


J2312-0109 


1.4262 


2006 


Oct 04 


4.3 


4.0 


J2340-0053 


1.3603 


2007 


Jan 21 


5.8 


3.9 






2007 


Apr 03 


5.8 


I qO.25,128 


J2358~1020 


1.1726 


2006 


May 23 


3.9 


3.6 






2006 


Jul 30 


5.5 


I gO.5,128 



Col. 1: Source name; col. 2: Mg II absorption redshift; col. 3: 
Date of observation; col. 4; Time on source; col. 5: Channel 
width in kms~^. Superscripts to channel width give the 
baseband bandwidth in MHz and number of spectra channels 
for the observation when a baseband set-up different from 
1 MHz split into 128 channels was used. 



ranta et al. 2005). We detect a weak 21-cm absorption to- 
wards this quasar with absorption redshift 2abs= 1.407, only 
~3000 km s~^away from the QSO emission redshift. Absorp- 
tion was confirmed by repeated observations at epochs sep- 
arated by seven months. We do not note any strong flux 
variation between these observations. The spectrum pre- 
sented in Fig. |4] results from the combination of all data. 
The absorption proflle is well fitted with a single Gaussian 
component with rp=0.008±0.001 and Aij=ll±2km s"\ 
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Figure 6. Apparent optical depth of 21-cm absorption features 
detected in our GMRT survey. The outermost dashed vertical 
lines give the velocity range over which the optical depth is inte- 
grated. The dotted vertical lines correspond to the width of the 
profile minus 5% of the absorption on both sides (see Fig. 1 of 
Ledoux et al. 2006). 



This source is unresolved at 8 GHz from MERLIN obser- 
vations part of the CLASS/JVAS survey with a beam size 
of 0.227" xO. 203". VLBA observations at 2 and 8 GHz re- 
solve the source into two components with a separation of 2 
and 1 mas respectively (Fey & Chariot 2000). The strongest 
component accounts for 85% of the flux detected in the 
VLBA map. However the total flux density 480 mjy mea- 
sured at 8 GHz in the VLBA observations accounts only for 
55% of the flux detected in the MERLIN image. It is not 
clear whether this difi^erence is due to flux variations or to 
a considerable fraction of the flux being emitted in a dif- 
fuse component that is resolved in the VLBA images. The 
H I column density derived from our GMRT observations is, 
iV(H l)=2.2xlO^^(Ts/100)(0.85//c)cm-2. Only Mg ll, Fe ll 
and Mg I absorptions are detected in the SDSS spectrum 
(See Fig. [71). 

J0850+5159: Radio source J0850-(-5159 is compact at 
mas scales in the VLBA images at 5 and 15 GHz (Tay- 
lor et al. 2005) implying /c -1.0. Strong (Tp=0.50) 21- 
cm absorption is detected towards this QSO which is ex- 
tremely red at optical and NIR wavelengths (u-K~4.8mag). 
The total 21-cm H I column density, with /c = 1, is 
iV(H l) = (2.81±0.09)xl0^''rsCm-2. Numerous absorption 
lines are detected in the SDSS spectrum in addition to Mg ll, 
Fe II and Mg I (see Fig. [7]). A 2175 A feature is also present at 
2.7(7 significance level at the redshift of the intervening Mg ll 
system (Srianand et al. 2008b) . An independent estimate of 
N{}i{ l) was derived by fitting a LMC2 extinction curve to 
the QSO spectrum, N{}1 1) = 5.73±1.10x 10^^ cm'^ There- 
fore we conclude that Ts ~ 204 ± 40 K (see Srianand et al. 
2008b for more details). 
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J0852+3435: Broad and shallow 21-cm absorption, con- 
firmed by observations at three different epochs, is detected 
towards this red QSO (u-K~5.6mag). Details about the 
2175 A feature detected at > 5a significance level at the 
redshift of the Mg ll absorption system are presented in 
Srianand et al. (2008b). Numerous absorption lines are de- 
tected in the SDSS spectrum in addition to Mg ll, Fe ll 
and Mg I (see Fig. [7]). The background radio source is com- 
pact in CLASS/JVAS and has flat radio spectrum (Jackson 
et al. 2007). Total 21-cm Hi column density with /c = 1 
is, iV(H l) = (1.27±0.06)xl0^^rscm"^. An independent esti- 
mate of A'^(H l) was derived by fitting the LMC2 extinction 
curve to the QSO spectrum iV(H l) = 6.97±1.30x 10^^ cm'^. 
Therefore we conclude that Ts ^ 550 ± 106 K (Srianand et 
al. 2008b) 

J1058-I-4939: The 21-cm absorption line detected towards 
this QSO is fitted well with two Gaussian components. The 
radio source is unresolved in the FIRST and GMRT images. 
The quasar is known to show variability at 151 MHz. Riley 
& Warner (1994) found that roughly 60% of the total fiux 
at 5 and 8 GHz comes from a region smaller than < 0.1" 
implying /c=0.6. VLBA or MERLIN images are not avail- 
able for this source. Integrated 21-cm optical depth yields 
Ar(Hi) = 1.25xlO^°(Ts/100)(0.60//c)cm-^ The SNRofthe 
SDSS spectrum is poor and only a few absorption lines are 
detected (see Fig. [7)l . 

J1623-f0718: 21-cm absorption is detected in two well de- 
tached components. The background radio source is unre- 
solved in GMRT and FIRST images. The source is com- 
pact even in CLASS/JVAS and has a fiat radio spec- 
trum (Jackson et al. 2007). Thus /c could be close to 1 
at ^100 mas scales. The integrated 21-cm optical depth 
yields Af(Hl) = 1.7xlO^**Ts cm'^ The Zn lH-Crll blend at 
A=2062 A, the Si ll and the Mn ll absorption lines are clearly 
detected even in the SDSS low dispersion spectrum (see 
Fig.0. 

J2340-0053: Radio source J2340-0053 is compact at mas 
scales in VLBA observations at 2 and 8 GHz (Kovalev et al. 
2007). The 21-cm absorption is narrow and the spectrum 
presented in Fig. |4] was obtained with a resolution of 1.0 
km s~^ per channel. The main component has a FWHM 
of 2.5±0.5 km s~^ which corresponds to an upper limit of 
136 K for the kinetic temperature. Clearly the gas responsi- 
ble for this system is cold. It is most likely that /c ~ 1 for this 
system. Under the assumption that /c = 1 and T < 200 K, 
we derive A'^(H l) < 9x10^^ cm~^. This means the compo- 
nent responsible for this 21-cm absorption is most probably 
a sub-DLA. An optical spectrum of this QSO was analysed 
by Khare et al. (2004). In Fig. |8]we show absorption lines 
associated to this systems on a velocity scale. The corre- 
sponding high spectral resolution KECK spectrum was ob- 
tained by Prochaska et al. (2007). The vertical dashed lines 
mark the location of the 21-cm absorption. It is clear that 
the 21-cm absorption is not originating from the component 
that dominates the equivalent width of strong metal lines. 



^ Note errors in Ts measurements for this system and the system 
a-t 2abs= 1.326 towards J0850-I-5159 are slightly reduced com- 
pared to that given in Srianand et al. (2008b) because of different 
procedures (we use here the errors obtained from the apparent 
optical depth profile technique). 




Relative velocity (km s"') 



Figure 8. Absorption lines in the 2abs= 1.3609 system towards 
J2340— 0053 observed at high spectral resolution (Prochaska et 
al. 2007). The dotted line spectrum corresponds to Al inA1862. 
Dashed vertical lines show the location of the 21-cm absorption 
line. 

Weak metal lines such as Fe llA2249 suggest that ~ 30% of 
the metals are in this component. The compactness of the 
background quasar and the fact that the metal absorption 
component associated with 21-cm is well detached from the 
central blend makes it an interesting system for measuring 
the time evolution of fundamental constants. However, the 
publicly available KECK data is not sufficient for this pur- 
pose. 

J2358— 1020: The background radio source is unresolved in 
high frequency VLBA observations (Fomalont et al. 2000; 
Fey & Chariot, 2000) and /c could be close to unity. 
The absorption profile is well fitted with a single Gaus- 
sian component whose FWHM implies a kinetic tempera- 
ture < 800 K. This gives log[iV(H l)]< 20.5 for = 1. The 
Zn ll-|-Cr 11A2O62A blend is clearly detected in the low dis- 
persion SDSS spectrum (Fig. 0. 



6 DETECTABILITY OF 21-CM ABSORPTION 

In this Section, we investigate the relationship between the 
detectability of 21-cm absorption and properties of UV ab- 
sorption lines from the Mg 11 systems. This may help in 
particular to preselect more efficiently candidates for future 
21-cm surveys. 

6.1 Rest equivalent width of Mg II absorption 

First we concentrate on Wr(Mg llA2796), the only quantity 
that has been used to define our present sample. In panel 
(a) of Fig. [9] we plot W^r(Mg llA2796) against the integrated 
21-cm optical depth for all the systems in our GMRT sam- 
ple. Detections (resp. non detections) are marked by filled 
(resp. open) circles. Even though the system with strongest 
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Figure 9. Top panels: The integrated 21-cm optical depth is plotted against various parameters of the UV absorption lines derived from 
the SDSS spectra. In the case of detection (filled circle) the optical depth is obtained by integrating over the observed absorption profile. 
Open circles are for 21-cm non-detections in which case limits are obtained by integrating the optical depth over a Gaussian component 
with peak optical depth corresponding to the 3a rms limit in the continuum and width 10 km s~^. Bottom panels: Different parameters 
of UV absorption lines are plotted against one another. Symbols are as in upper panels. The size of the symbols are scaled by the 21-cm 
optical depth limits (per lOkms"^) reached. The vertical dashed line in all the panels gives the median value. Other lines indicate the 
limiting values and/or allowed ranges as discussed in the text. 



21-cm absorption in our sample has the largest equivalent 
width, we do not find evidence for a correlation between 
the two quantities. Indeed, the median Wr(Mg nA2796) is 
1.85 A for the systems observed with GMRT and there are 
only three (out of 9) 21-cm detections from systems having 
H/r(MgiiA2796)>1.85A. As we have reached similar limiting 
21-cm optical depth over the whole range of Wr(Mg nA2796) 
sampled by our survey, this can not be due to any artifact 
introduced by variation of the optical depth limits through 
the sample. 

For comparison with our sample, we summarize in Ta- 
ble [S] the properties of systems with Wr(Mg nA2796)> 
1 A in the low-z sample of Lane (2000). The median 
Wr(Mg nA2796) and redshift in this subsample is 1.75 A 
and 2abs=0.7. There are 2 and 5 detections respectively in 
systems with Wr(Mg nA2796) below and above the median 
value. It seems therefore that the probability to detect 21- 
cm absorption in strong Mg ll systems is higher at lower 
redshift. We further investigate this issue in the next Sec- 
tion. 



6.2 Mg 11 Doublet ratio (DR) 

In panel (b) of Fig. |9] we plot the integrated 21-cm optical 
depth vs. the Mg ll doublet ratio (DR). It is clear that 21- 
cm detected systems are confined to regions with DR<1.3. 
For the systems with Wr(Mg nA2796)> lA and DR<1.3 
we detect 21-cm absorption in roughly 33^11% (8/24) of 
the systems to be compared to the overall success rate in 
the sample, 24+g^% (8/33). The median value of DR m our 



sample is 1.12 and only two systems with DR greater than 
this are detected in 21-cm absorption. 

As pointed out by Petitjean & Bergeron (1990) the 
number of Voigt profile components found in high resolu- 
tion spectra correlates well with the total equivalent width 
of blended lines (see Fig. 3 of Ellison 2006). Thus equivalent 
widths measured in SDSS data reflect primarily the num- 
ber of individual components and/or the velocity spread 
between them and not directly the column density (see 
P06). On the contrary, a doublet ratio close to 1, even in 
a low dispersion spectrum, means high level of saturation 
and therefore high column density. Thus systems with high 
H^r(Mg nA2796) and DR~1 provide higher probability of a 
line of sight hitting high Mg ll column density clouds. 

Interestingly, in the GMRT sample, there are 9 systems 
with DR>1.3 and 7 of these systems have W^r(Mg nA2796) 
> 1.85 A (see panel e of Fig. |9]). The lower detection rate 
of 21-cm absorption in systems with Wr(Mg nA2796)> 1.85 
in the GMRT sample could be due to high values of DR in 
most of these systems. It is most likely that the high value of 
VFr(Mg ll) in these systems is primarily due to the velocity 
field and not to a large Mg ll column density. 

From Fig. 3.3 of Lane (2000) we find that all the low- 
z 21-cm detections have DR<1.4. Thus our finding of high 
detection rate of 21-cm absorption in systems with low DR 
is consistent with what is seen for low-z 21-cm absorbers. 
The difference between high and low redshift sample noted 
in the previous Section is therefore probably related to a 
redshift evolution in the overall structure of Mg ll systems. 
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6.3 The Wr(Mg llA2796)/Wr(Fe llA2600) ratio 

RTN06 showed that DLAs have 1.0< 
H/r(Mg iiA2796)/Wr(Fe iiA2600)< 2. We shall caU this 
ratio _Ri. In their Mg ll sample there are 19 systems in the 
redshift range 1.10 < z < 1.45, with Wr(Mg nA2796)> lA 
and Ri in the above mentioned range. Seven of them (i.e 
37ti^%) have log N{Y{ l)>20.3 and are DLAs as per the 
conventional definition. However, 18 of these 19 systems 
have log 7V(H l)>19.6. 

From Panel (c) of Fig|5]it is clear that our selection of 
systems based on W^r(Mg iiA2796) alone has ensured that a 
large fraction of systems (24/33) in our sample have _Ri in 
the above mentioned range. In addition, the detection rate 
of 21-cm absorption in our sample is high when Ri is low. 
The median value of R\ in our sample is 1.54. About 411^5% 
(resp. llt'fVo) of the systems with Ri < 1.54 (resp. >1.54) 
show 21-cm absorption. In the -Ri range covered by DLAs 
in RTN06 we find that there are 24 Mg ll systems in our 
sample with 7 of them (i.e 29tii%) showing detectable 21- 
cm absorption. The other two detections have Ri consistent 
with 2 within measurement uncertainties. For comparison, 
in the low- 2; sample of Lane (2000) there are 11 systems with 
1 < _Ri < 2 and 6 of them show 21-cm absorption. 

In panel (f) of Fig. |9]we plot Ri vs DR. Out of the 20 
systems with both DR<1.3 and i?i < 2, 7 systems show 21- 
cm absorption (i.e 35li[3%) systems. This is very close to the 
detection rate we obtained with the condition DR<1.3 only. 
Reducing the upper limit of Ri to 1.6 produces an increase 
of the detection rate to 41^jj%. 

Finally we notice that roughly half the systems with 
Wr(Mg nA2796)> 1.85A also have Ri > 2. This again sug- 
gests that in this redshift range large Wr(Mg llA2796) does 
not guarantee high column density. 

6.4 The Wr(Mg iA2852)/VKr(Mg nA2796) ratio 

From the data listed in Table 1 of RTN06 we find 
that the DLA population with H^r(Mg nA2796)> lA has 
Wr(Mg iA2852)/iy"r(Mg nA2798) (called R2 from now on) 
in the range 0.01 to 0.6. Conversely, ^30^6% of systems 
with 0.01 < R2 < 0.6 are DLAs. The fraction increases to 
^52tll% for R2 > 0.3. In addition, 92t%% of the systems 
with R2 > 0.3 have log 7V(II l)> 19.60. From the sample of 
Lane (2000) we find that 35% (resp. 40%) of the systems 
with Wr(Mg llA2796)>lA and R2 > 0.1 (resp. R2 > 0.3) 
show detectable 21-cm absorption. Note however that the 
Mg lA2852 equivalent width is not always available. 

It is clear from panel (d) of Fig. [9] that all our 21-cm 
detections, in consistency with the results from low-redshift 
studies, are associated with R2 > 0.1. We find that the de- 
tection rate of 21-cm absorption is ~36ti2% when R2 > 0.1. 
About 44+jg% of the systems with R2 >0.2 (the median 
value in our sample) show 21-cm absorption. If we restrict 
ourself to R2 > 0.3 the detection rate in our sample increases 
to 70^26^- It is also apparent from panel (g) of Fig [9] that 
a joint constraint, DR< 1.3 and R2 > 0.3 (solid box in this 
panel), gives a 88^31% detection rate (7 out of 8) for 21-cm 
absorption. Thus it appears that using joint constraints on 
DR and R2 can yield a very high detection rate of 21-cm ab- 
sorption. This is confirmed by the low-z data of Lane(2000) 
summarized in Table |6] There are only 5 systems (exclud- 



Table 5. Sample of sources with milli-arcsecond scale images 



QSO 
(1) 


fc 

(VLBA) 
(2) 


LAS 
(mas) 
(3) 


VKr(Mg IIA2796) 

(A) 

(4) 


1 r I 
-J ^"'^ 
(km s-i) 

(5) 


J0108-0037 


0.30-0.701 


5.5 


0.43±0.05 


1.84-4.3 


J0240-2309 


0.65^ 


10.0 


1.85±0.002 


< 0.03 


J0259-0019 


i.ooi 


<15 


1.76±0.04 


< 0.28 


J0742-I-3944 


0.58^ 


8.1 


2.48±0.41 


< 0.62 


J0748-I-3006 


0.413 


54.3 


3.68±0.D6 


< 0.44 


J0808+4950 


0.85^ 


9.7 


1.33±0.09 


0.12 


J0850-I-5159 


0.993 


3.0 


4.56±0.12 


15.5 


J0953-I-3225 


0.843 


3.2 


1.55±0.05 


< 0.26 


J1126-I-4516 


0.573 


5.6 


1.26±0.25 


< 0.14 


J1208-I-5441 


0.883 


1.2 


1.85±0.25 


< 0.13 


J1508-I-3347 


0.743 


11.4 


2.67±0.08 


< 0.31 


J2340-0053 


0.87-* 


13.6 


1.60±0.04 


2.87 


J2358-1020 


0.98^ 


2.2 


1.52±0.15 


0.24 



1 Beasley et al. 2002 (2 & 8 GHz); ^ Fomalont et al. 2000 (5 
GHz); 3 Helmboldt et al. 2007 (5 GHz); " Kovalev et al. 2007 (2 
& 8 GHz) 



ing poor upper limits on Mg l) with R2 > 0.3 and 3 of 
them show detectable 21-cm absorption. One of these non- 
detections (zabs= 0.9004 towards 1629-1-120) has DR = 1.6 
and the other system (zabs^ 0.7713 towards 1556—245) has 
Ri > 5. 

From panel (h) of Fig. |9] it can be seen that joint con- 
straints on Ri and R2 does not produce higher detection 
rate compared to that obtained from the joint constraints 
on DR and Ri. 



6.5 MilHarcsecond scale radio structure 

We address here possible effects introduced by issues related 
to covering factor of the background source that may play 
an important role on the detectability of 21-cm absorption 
Although VLBI spectroscopy at 467 MHz of the Zabs= 2.04 
DLA towards PKS 0458-020 suggests that the absorbing gas 
covers a large fraction of the extended (>17kpc) radio com- 
ponents (Briggs et al. 1989), time- variability of 21-cm op- 
tical depth seen in the absorbers towards B0235-)-164 and 
B1127— 145 suggest that absorbing gas is patchy at sub-kpc 
scales (Wolfe, Briggs & Davis, 1982; Kanekar & Chengalur, 
2001). In addition, spatially resolved spectroscopy of the 
three closely spaced (<1 arcsec) images of APM 0827-1-5255 
have revealed a strong variation in Mg 11 optical depths over 
kpc scales (Petitjean et al. 2000; Ellison et al. 2004). 

At present, low frequency (i.e <1.4 GHz ) mas scale 
images are not available for any of the sources in our sam- 
ple. Our ongoing VLBA observations at low frequencies will 
provide better constraints on fc in the near future. However, 
high frequency (i.e > 2 GHz) milli-arcsecond scale images 
are available for 13 sources out of the 22 sources that are 
compact in FIRST and our 610 MHz images. We give details 
of these objects in Tabled Measured covering fraction (/c) 
and largest angular size (LAS) are given in columns 2 and 3 
respectively. Covering fraction has been estimated as the ra- 
tio of the core fiux density to the total flux density detected 
in the VLBA images. To achieve as much uniformity as pos- 
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Table 6. Sample of Mg n systems with Wr(Mg llA2796)>lA from Lane (2000). 



Source name 






Wr(Mg IIA2796) 


Wr(Mg IIA2803) 


Wr{Mg IA2852) 


Wr(Fe IIA2600) 


dv 




/rdv 








(A) 


(A) 


(A) 


(A) 


(kms^-"-) 




(kms-l) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


0109+176 


2.155 


0.8392 


1.75 


1.20 


<0.20 


1.09 


9.1 


0.108 


<1.09 


0109+200 


0.746 


0.5346 


2.26 


1.71 






7.6 


0.029 


<0.27 


0229+341 


1.240 


0.7754 


1.92 


2.02 


<1.13 




8.8 


0.043 


<0.43 


0235+164 


0.940 


0.5238 


2.42 


2.34 


0.91 


1.79 


1.0 + 


0.007 


13.0 


0248+430 


1.316 


0.3939 


1.86 


1.42 


0.70 


1.03 


5.8 


0.032 


2.99 


0420—014 


0.915 


0.6330 


1.02 


0.86 


<0.36 


.... 


10.6 


0.146 


<1.59 


0454+039 


1.343 


0.8597 


1.53 


1.40 


0.37 


1.11 


10.0 


0.010 


<0.11 


0735+178' 




0.4246 


1.33 


1.04 


0.18 


0.87 


3.5 


0.058 


<0.36 


0827+243 


0.939 


0.524 


2.90 


2.20 




1.90 


10.0* 


0.004 


0.22 


0805+046 


2.876 


0.9598 


1.01 


0.83 


<0.60 


0.24 


4.9 


0.095 


<0.71 


0957+003 


0.907 


0.6720 


1.77 


1.32 






8.3 


0.105 


<1.01 


1127-145 


1.187 


0.3130 


2.21 


1.90 


1.14 


1.14 


5.4 


0.021 


3.11 


1218+339 


1.519 


0.7423 


1.34 


1.08 


<0.70 


1.00 


8.6 


0.020 


<0.20 


1229-021 


1.038 


0.3950 


2.22 


1.93 


0.49 


1.59 


i.ot 


0.007 


3.00 


1327-206 


1.169 


0.8530 


2.11 


1.48 


<0.40 


0.76 


9.2 


0.201 


<2.04 


1354+258 


2.004 


0.8585 


1.00 


0.86 


<0.10 


<0.20 


9.2 


0.103 


<1.05 


1556-245 


2.815 


0.7713 


2.07 


1.91 


1.07 


<0.20 


8.8 


0.068 


<0.68 


1622+238 


0.927 


0.6560 


1.29 


<1.69 


<0.40 


1.13 


10.0* 


0.004 


2.20 


1629+120 


1.795 


0.5313 


1.40 


1.35 


0.31 


0.70 


0.6* 


0.004 


0.49 






0.9004 


1.06 


0.67 


0.44 


0.63 


9.4 


0.080 


<0.82 


2212-299 


2.703 


0.6329 


1.26 


1.00 


0.36 




8.1 


0.217 


<2.07 



Col. 1: source name as given in Lane (2000); col. 2: QSO emission redshift; col. 3: Absorption redshift of Mg II system; cols. 4, 5, 6 and 
7: Rest equivalent widths of Mg llA2796, Mg llA2803, Mg lA2852and Fe llA2600; cols. 8 and 9: spectral resolution and 3a uncertainty on 
the optical depth; col. 10: integrated 21-cm optical depth for the detections or 3a upper limit to it in case of non-detections. The upper 
limits on J rdv have been computed for dv=10kms^^. 

t BL Lacertae object a featureless continuum except for the intervening absorption line system at Zijbs=0-4246 (Carswell et al. 1974). 
t 21-cm absorption data for 0235+164, 1229-021 and 1622+238 are originally from Wolfe et al. (1978), Brown & Spencer (1979) and 
Curran et al. (2007) respectively. 

* 21-cm absorption data for 0827+243 and 1629+120 are from Kanekar et al. (2003). 




Figure 10. The ratio Wr(Mg lA2852)/iyr (Mg IIA2796) vs. fc. 
The filled and open circles are for systems with and without 21- 
cm absorption respectively. 



sible, we have used the VLBA 5 GHz (rather than VCS 2 
and 8 GHz) images whenever available. From this table it 
appears that the detectabiUty of 21-cm does not depend on 
the covering factor in this sub-sample. We plot in Fig. [10] 
the ratio R2 against the measured covering factor fc for sys- 
tems with (filled circles) and without (open circles) 21-cm 
detection. It is apparent that while the distribution of /c is 
similar between detections and non-detections, R2 separates 
the two very nicely. 

A covering factor less than 1 increases the Scr optical 
depth limit by a factor of In order to quantify the de- 

tection rate of 21-cm absorption we give in the last column 
of Table [5] the integrated optical depth corrected from par- 
tial coverage. We come back to this subsample in the next 
Section while discussing the redshift distribution of 21-cm 
absorbers. 



7 REDSHIFT DISTRIBUTION OF 21-CM 
ABSORBERS 

The number density per unit redshift 7121 (To, Wo, z), of 21- 
cm systems with integrated optical depth T21 > To and Mg 11 
equivalent width Wr > Wo can be obtained from the number 
density per unit redshift of Mg 11 absorbers, nMgii(Wo,2), 
using the equation. 
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Table 7. Rcdshift distribution of 21-cm absorbers 



0.3 



Sample 




Wo 


% 


N 


N21 


C 




"21 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 




(7) 


GMRT 




1.0 


0.10 


9 


4 


44+0-35 





4+0.094 
-0.057 


(-^abs 


1.3) 




0.30 


26 


3 


12+0" 





r,on+o.030 

"■^"-0.017 








0.50 


33 


4 







nqi +0.026 

"•"■^-o.oie 








0.50" 


12 


2 


Q ^7+0.22 
'-'•^'-0.11 





" -0.028 






1.8 


0.30 


14 


1 


'-'•'^'-0.06 





006+0-013 
'^^°-0.005 








0.50 


20 


2 


10+0" 





nns+oolo 

UU«_0.006 


Low- 2 




1.0 


0.30 


9 


5 







n7Q+0.055 
""'-0.035 


(-^abs 


0.5) 




0.50 


11 


4 


-^6+0-29 





051+0-042 
'^^1-0.025 






1.8 


0.30 


6 


4 


n R7+0'53 

'-'•°'-0.32 





r,r,o+o.020 

'J^'^-0.012 








0.50 


7 


4 


n 57+0 « 





^^'-'-0.011 



Col. 1: Sample name and mean redshift; col. 2; Equivalent width 

limit; col. 3: Threshold of J Tdv in kms~i; col. 4: Number of 

systems in the sample with J T3a-_iodv< To; col. 5: Number of 

21-cm absorption detections with integration line depth 

J rdv> To; col. 6: Fraction of Mg II systems with detectable 

21-cm absorption; col. 7: Number density per unit redshift of 

21-cm absorbers with J Tdv> To. 

" considering only the subsample in Table [5] 



n2i(T2i > To, Wr >Wo,z)^Cx nMgii(Vyr > Wo, z). (3) 

Where, C is the fraction of Mg 11 systems with 
Wr{Mg llA2796)>'Pyo that show detectable 21-cm absorp- 
tion. 

The redshift incidence of Mg 11 absorbers can be ap- 
proximated as, 



nMgii(W^o, 2) = no X (1 + zy 



(4) 



with no = O.O8OI0.005 and 7 = 1.40 ± 0.16 for Wo = 
lA and no = O.Oie+o.ws and 7 = 1.92 ± 0.30 for Wr > 
I.SA (Prochter et al. 2006). Therefore at <2abs> = 1.3, 
the mean redshift of our sample, we derive, nMgii(Wo = 
lA,z = 1.3) = 0.257toml and nMgii(H/o = 1.8A,z = 1.3) 

- 079+0-054 

— u.u(y_oo3o- 

The fraction of Mg 11 systems that show detectable 21- 
cm absorption with integrated optical depth, T21 > To, 
is estimated as follows. First we consider only the sys- 
tems (say N) with radio spectra sensitive enough so that 
J r3CTdv< To- Since the typical resolution of the Lane's spec- 
tra is 7-10 km s~^ (see below), we have computed the values 
of the 21-cm 3cr-detection limit for a velocity resolution of 
lOkms"^ (see Table [3]). Then the number of these systems 
with detectable 21-cm absorption (referred to as N21) is es- 
timated by considering the detections that have integrated 
line depth, 72i > Tq. Then C is given by N21/N and the er- 
ror in C is computed using small number Poisson statistics 
(Gehrels, 1986) . The value of n2i for various sub-samples 
defined using To and Wo are summarised in Table [T] For 
Wo = 1 A the sub-sample with To = 0.5 contains 33 Mg 11 
systems. About 12% of these systems are 21-cm absorbers. 
If we restrict ourselves to compact QSOs with VLBA mea- 
surements, and use To corrected for partial coverage (see 
Table. [Sj, we find that 17% of these systems are 21-cm ab- 
sorbers. Thus uncertainties due to unknown covering factor 
should have little effect on our conclusions. 




Figure 11. Number of 21-cm absorbers per unit rcdshift, 
?i2i(z), for integrated 21-cm optical depth, T > 0.3 and 
Wo(Mg 11) > I.OA (solid symbols) and VKo(Mg 11) > 1.8 A (open 
symbols). Lines and hashed areas show the number of absorbers 
per unit redshift for DLAs (RTN06, dotted line) and Mg II ab- 
sorbers with Wo >lA (P06, dashed line). 



We use the sample of Lane (2000, see Table |6]l to study 
the redshift evolution of n2i . For two of the systems in this 
sample (zabs^ 0.524 towards 0827-f 243 and Zabs^ 0.5313 to- 
wards 1629-1-120) we use the results of better quality GMRT 
observations by Kanekar et al. (2003). In the case of the sys- 
tem at 2abs= 0.6560 towards 1622-1-238, the 21-cm absorber 
detected by Curran et al. (2007) is most probably located in 
front of the extended lobes. As we have avoided such systems 
in our sample, we do not include it in our statistical analy- 
sis. The values of n2i for this sub-sample with < z >~0.5 
are given in Table [7] 

From Table 4 of Prochter et al. (2006) we expect the 
number density of Mg 11 systems with Wo ~ 1.0 and 1.8 A to 
increase by a factor of 1.8 and 2.3 respectively between z — 
0.5 and 1.3. If physical conditions in Mg 11 absorbers (like 
A'^(H l) and CNM fraction) do not change with redshift we 
expect C to remain constant and therefore n2i to increase 
with redshift. This is not the case however. It is clear from 
Table [7] that, for a given value of Wo and To, C is less at 
higher z. This decrease is even more pronounced for higher 
equivalent width (VKr(Mg llA2796)>1.8A). As the decrease 
in C is stronger than the increase in the number of Mg 11 
systems per unit redshift interval we find n2i decreases with 
redshift. 

Rao, Turnshek & Nestor (2006) give the number of 
DLAs : nDLA(2 = 0.609) = 0.079 ± 0.019 and nm,A{z = 
1.219) = 0.120 ± 0.025 suggesting a factor 1.5 increase in 
the number per unit redshift of DLAs. The number of 21- 
cm absorbers with Wo > lA and To > 0.3 matches well the 
number of DLAs at z ~ 0.5 but it is a factor of 4 smaller at 
z - 1.3 (see Fig. [Til- 
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Table 8. Results of non-parametric and rank correlation tests. 



2 3 
W^(MgII \2796) 

Figure 12. Width of the 21-cm absorption as a function of the 
Mg IIA2796 rest equivalent width. Filled hexagons are for our 
GMRT sample, open triangles for the low-z: sample and open 
hexagons for z > 1.8 DLAs. 



The decrease in n2i with redshift can be a consequence 
of the evolution of H I column density in the Mg ll absorbers 
or lower CNM covering factor at high redshift which may re- 
flect a strong evolution either in the filling factor of the neu- 
tral cold gas or the fc of the radio sources. From Table [T] it 
is apparent that the fraction of 21-cm absorbers is not larger 
for Wr{Mg iiA2796)>1.8 A compared to Wr{Mg iiA2796)>l 
A at z ~ 1.3. This is counter-intuitive as RTN06 have 
found that the fraction of DLAs amongst Mg ll systems in- 
creases with increasing Wl (Mg ll) (see their Fig. 6). Now the 
% = 0.3 corresponds to log Af(H l) = 19.74 for the Ts = 100 
K. From the sample of RTN06 we find that the fraction of 
Wr > lA systems with log iV(H l) > 19.74 for the 2abs<1.0 
and Zaba>1.0 is 50±10% and 75±15% respectively. Thus the 
decrease in n2i for Ta — 0.3 is unlikely to be due to the evolu- 
tion in column density amongst the Mg ll absorbers. Further 
we have shown that the fraction of 21-cm absorbers in our 
sample does not depend strongly on radio morphology (see 
Table [7]) • The available data are therefore consistent with 
the CNM fraction of neutral H I in Mg ll systems being 
smaller at higher redshifts. This will also be consistent with 
the results of Menard et al. (2008) who found the dust con- 
tent of the strong Mg ll absorbers to be increasing with the 
cosmic time over the redshift range of 0.4 < z < 2.2. 



8 VELOCITY SPREAD OF 21-CM 
ABSORPTION 

In this Section we investigate any possible correlation be- 
tween the velocity width of the 21-cm absorption and the 
Mg IIA2796 rest equivalent width. Indeed, Lane (2000) sug- 
gested that systems with large metal line equivalent widths 



Sample* 


P 


T 


Probability 




(1) 


(2) 


(3) 


(4) 


(5) 




GMRT(9) 


0.133 


0.056 


0.732 


0.835 


0.4 


GMRT-J0804+3012(8) 


0.429 


0.214 


0.289 


0.457 


1.0 


GMRT -|-high-z(12) 


0.165 


0.092 


0.608 


0.676 


0.6 


low-2(10) 


0.596 


0.494 


0.069 


0.046 


1.8 


low-2 -0809-1-483 


0.678 


0.592 


0.045 


0.026 


2.0 


All(22) 


0.482 


0.337 


0.023 


0.028 


2.3 


All-0804-0809(20) 


0.573 


0.411 


0.008 


0.011 


2.8 



* Numbers in brackets correspond to the number of 21-cm 
absorbers in the subsamples. 



could be associated with multi-component 21-cm absorp- 
tion. In addition, using mostly the sample described in Lane 
(2000), Curran et al. (2007) have reported a correlation (sig- 
nificant at the 2.2—2.8(7 level) between the 21-cm line width 
and the Mg llA2796 equivalent width. 

In order to probe the presence of such a correlation in 
our sample we computed the 21-cm velocity width using the 
apparent optical depth profile and the method described in 
Ledoux et al. (2006) for the low ionization metal lines (see 
Section 5 and Fig. [6]). Even though our parent sample of 
Mg II systems has a wide range in Wi (Mg llA2796) most of 
our 21-cm detections correspond to systems with widths of 
1-2 A. However, we find that measured 21-cm absorption 
velocity width ranges from 7 to 130 km/s (Table 4). 

In Fig. [12] we plot the 21-cm line width as a function of 
Wr{Mg IIA2796). Filled hexagons are for our GMRT sample, 
open triangles for the low-z sample and open hexagons for 
z > 1.8 DLAs. Results of non-parametric and rank correla- 
tion tests are summarized in Table [8] p and r are the Spear- 
man and Kendall rank correlation coefficients respectively. 
The probability for the measured correlation to happen by 
chance is given in the 4th and 5th columns. The last col- 
umn gives the significance of the correlation assuming the 
statistics to be Gaussian for the spearman rank correlation 
test. 

We do not find any significant correlation between Wi 
and the 21-cm velocity width in the GMRT sample. The 
largest 21-cm velocity spread is seen towards J0804-I-3012. 
Since the velocity spread could be due to the complex struc- 
ture of the radio source that is resolved in our GMRT map 
(see discussions in Section 5) we removed this system from 
the analysis. The possible correlation is stronger but still 
not significant (^la assuming statistics to be Gaussian). 
Adding 3 high-z 21-cm absorbers (open hexagons in Fie:. ll2|l 
observed at 1.8 < z < 2.1 for which Mg ll equivalent widths 
are known does not modify the original conclusion of no 
correlation (^Ict). These 3 21-cm absorbers are Zabs=2.0394 
towards 0458-020 (Wolfe et al. 1985; Wolfe et al. 1993), 
Zabs=1.9440 towards 1157-^014 (Wolfe & Davis 1979) and 
Zabs=1.7763 towards 1331 + 170 (Wolfe & Davis 1979; Briggs 
et al. 1983). 

The results for the low-z data are plotted as open tri- 
angles. Here, we use the 6 21-cm absorbers from the low-z 
sample discussed in the previous Section (see Table. [6]) plus 
4 additional 21-cm absorption systems with Mg ll equiva- 
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lent width measurements. These four additional systems are 
«^bs=0.2213 towards 0738+313 (Lane et al. 1998; Kanekar 
et al. 2001), 2abB=0.4368 towards 0809+483 (Aldcroft et al. 
1994; Briggs et al. 2001), Zabs=0.2377 towards 0952+179 
(Lane 2000; Kanekar ct al. 2001) and Zabs=0.692 towards 
1328+307 (Davis & May 1978; Cohen et al. 1994). The cor- 
relation in this low-2 sample (with or without 0809+483 
which has a 21-cm velocity width of 150kms~^) is ~2cr sig- 
nificant. There is therefore a difference between low and high 
redshift. The weak correlation (2—2.5 a level for a Gaussian 
distribution) noted by Curran et al (2007) is present when 
we consider the full-sample i.e. all the sources at low and 
high redshift together. 

The main differences between low- 2: and GMRT sam- 
ples are (i) a paucity of 21-cm absorption systems with 
Wr{Mg iiA2796)<l A and (ii) a lack of 21-cm absorbers 
among Mg II systems with Wr(Mg 11) > 1.8 A in the GMRT 
sample. While the first difference is probably related to our 
sample selection the second difi'ercncc is probably due to a 
true evolution in the physical state of the Mg 11 systems with 
redshift (see discussion in the previous Section). Therefore 
we think it is not wise to mix the absorbers at low and high 
redshifts. 

The absence of strong correlation between the total 
Mg II equivalent width and the 21-cm absorption width is 
consistent with the idea that the Mg 11 equivalent width is 
mostly correlated with the overall kinematics of the gas in 
the absorbing system and not with the column density in 
the component associated with the cold gas. In the case of 
J2340-0053 where we have access to high resolution echelle 
spectra, we find that the 21-cm absorption is not associ- 
ated with the strongest metal line component (see Fig. 9). 
In addition as pointed out by Gupta ct al. (2007) wc usually 
find that the 21-cm absorption redshift is slightly different 
(within 70 km s~^) from the mean redshift defined by the 
metal absorption linos in the SDSS spectrum. This can again 
be explained if the 21-cm absorption does not originate from 
the strongest metal line components that define the absorp- 
tion redshift at low resolution. 



9 SUMMARY AND DISCUSSION 

We have reported the results of a systematic GMRT sur- 
vey of 21-cm absorption in a representative sample of strong 
Mg II systems selected from the SDSS spectra in the redshift 
range 1.10< z <1.45. We have performed GMRT observa- 
tions of ~70% (resp. 35%) of our initial sample of Mg 11 
systems detected in DR5 with Wr{Mg llA2796)>l A located 
in front of compact radio sources with flux density > 100 mjy 
(resp. in the range 50-100 mJy). Wc show that our observed 
sample is representative of the full parent population. 

We report detection of 9 new 21-cm absorption systems 
out of 35 systems observed in our survey. This is by far 
the largest number of systems detected in a single system- 
atic survey in a narrow redshift range. Two of those sys- 
tems also show 2175 A dust feature at the redshift of the 
absorbers (Srianand et al. 2008b). We provide spectra of all 
non-detections with upper limits on the 21-cm optical depth. 

We study the dependence of detectability of 21-cm ab- 
sorption on different properties of the UV absorption lines 
detected in the SDSS spectra. We find that if absorption 



systems are selected with a Mg 11 doublet ratio, DR < 1.3, 
and a ratio Wr{Mg l)/Wr(Mg 11) > 0.3, the success rate for 
21-cm detection is very high (up to 90%). We notice that 
the detections found in a low-2 sample by Lane (2000) also 
obey these joint constraints. Low value of DR and high val- 
ues of Wr(Mg i)/IVr(Mg 11) in a system will mean a high 
column density of Mg 11 and high neutral fraction. Thus 
this selection favors objects with high N(H i). In our sam- 
ple, we find an apparent paucity of 21-cm absorption among 
systems with Wr(Mg llA2796)>1.8 A, the median Wi of our 
sample. This is contrary to what has been seen at \ow-z 
(Lane 2000). Interestingly most of these high Wt systems 
have high DR and low values of Wr{Mg l)/T-Vr(Mg 11). This 
strongly suggests that the equivalent width in these systems 
is dominated by velocity spread and not by line saturation. 

We discuss the number of 21-cm absorption systems per 
unit redshift interval for a given limiting value of the inte- 
grated 21-cm optical depth and Wr(Mg llA2796). Wc show 
that the fraction of Mg 11 systems with 21-cm absorption 
and the 7121 decrease from z ~ 0.5 to z ~ 1.3. The de- 
crease is larger when we use higher equivalent width cut- 
off. Using a sub-sample of compact sources, with high fre- 
quency VLBA observations available, wo show that this can 
not be accounted for by simple covering factor effects. As 
mentioned above and based on the available data, it ap- 
pears that most likely the main reason behind this cosmo- 
logical evolution is the decrease of the CNM covering factor 
(and volume filling factor) in the strong Mg 11 absorbers. 
Indeed, it is known that the number of Mg 11 systems per 
unit redshift increases with increasing redshift. The evolu- 
tion is steeper for stronger systems (Steidol & Sargent, 1992 
and P06 for recent reference). Using the data of Steidol & 
Sargent (1992), Sriamand (1996) found that the strongest 
redshift evolution was seen among the Mg 11 absorbers with 
Wr(Fe ii)A2383/Wr(Mg ii)A2796<0.5. This clearly means 
the physical conditions in strong Mg II absorbers are dif- 
ferent at high and low-z. 

Previous surveys (listed in Section 2) have hinted the 
lack of 21-cm absorption systems at high-2 {z > 2) compared 
to what is seen at low-2, 2 <0.5. The low detection rate of 
21-cm absorption in high-2 DLAs is attributed to either the 
gas being warm (high Tg) or to low values of the covering 
factor because of high-2 geometric effects (see Kanekar & 
Chengalur, 2003 and Curran & Webb, 2006, respectively). 
In two of our systems (at 2abs= 1.326 towards J0850+5159 
and 2abs= 1.209 towards J0852+3435), the presence of a 
2175A dust feature allows us to derive an estimate of A'^(H i) 
and therefore Tg (Srianand et al. 2008b). Both systems have 
high metallicity, large iV(H l) (~6xl0^^ cm"^) and % con- 
sistent with the temperature in cold neutral (CNM) gas. In 
the case of the system at 2abs= 1.361 towards J2340— 0053, 
the 21-cm line width is consistent with the absorbing gas 
having T < 200 K. Therefore, we already have 3 systems 
in our sample with properties consistent with the condi- 
tions expected for CNM gas. Only 1 such system has yet 
been detected at z > 2 (York et al. 2007). Thus our results 
are consistent with the increase of the CNM covering factor 
with decreasing redshift. However, to ascertain this result 
it would be important to derive A'^(II l) directly from the 
Lyman-Q absorption. Thus followup HST/COS observations 
of a UV bright subsample of our sample is very important 
to study the redshift evolution of Tg. As pointed out by nu- 
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merous authors, the physical state of the gas depends on 
the star formation activity in the vicinity of the absorber. 
Models of multiphase interstellar medium suggest that the 
temperature of different phases depends on the background 
UV radiation field, metallicity, dust content and cosmic ray 
density and therefore on the local star formation activity. 
It is possible to conjecture that the evolution of 7221 as a 
function of redshift is probably related to the fact that at 
z ~ 1.3 the average SFR in a galaxy is larger than that at 
z ~ 0.5 so that we expect the CNM fraction to increase with 
decreasing redshift. 

We selected systems with Wr > 1 A but detected by 
chance a 21-cm absorption in a system with Wr = 0.43 A(at 
2abs = 1.3710 toward J0108-0037). Efforts are underway at 
GMRT to extend our survey to weaker (Wr <lA ) Mg 11 
systems. This will be crucial for understanding the physical 
state of Mg 11 systems and to determine the detectability 
of 21-cm absorption versus Wr- Ideally one would like to 
estimate the number density of 21-cm absorbers and mea- 
sure the cosmological evolution without preselection from 
the UV absorption lines. This can be achieved only by a 
bhnd survey of 21-cm absorption in front of radio loud 
QSOs. It will be possible to embark upon such a survey 
with the upcoming Square Kilometer Array (SKA) pathfind- 
ers. In particular, the Australian Square Kilometer Array 
Pathfinder (ASKAP) with its instantaneous wide bandwidth 
of 300 MHz and large field of view (30 degree^) is an ideal 
instrument for this (Johnston et al. 2008). An ASKAP sur- 
vey with 150 pointings of 16 hrs each (i.e 2400 hrs in total) 
in the 700-1000 MHz frequency band would yield detection 
of ~100 to 250 intervening 21-cm absorbers in the redshift 
range 0.4< z <1. 

We have estimated the velocity spread of the 21-cm ab- 
sorptions using the apparent optical depth method (Ledoux 
et al. 2006). We do not find any statistically significant cor- 
relation between Wr(Mg llA2796) and the 21-cm velocity 
width in our sample. A marginal correlation is found for 
the low-z sample. The absence of correlation in the high- 
z sample is related to the lack of 21-cm absorbers among 
Mg II systems with Wr{Mg ll)>1.8 A in the GMRT sam- 
ple. This is probably due to a true evolution with redshift of 
the physical state of the Mg 11 systems and consistent with 
the idea that the Mg 11 equivalent width is mostly corre- 
lated with the overall kinematics of the gas in the absorbing 
system and not with the column density in the component 
associated with the cold gas. When high spectral resolution 
data are available, we note that the 21-cm absorption is not 
always associated with the strongest Mg 11 component. 

As the energy of the 21-cm transition is proportional to 
X = a^Gp/^, high resolution optical and 21-cm spectra can 
be used together to probe the combined cosmological vari- 
ation of these constants (Tubbs & Wolfe, 1980). Using this 
technique on scanned data, Tzanavaris et al. (2005) obtained 
Ax/x = (0.63 ± 0.99) X 10"^ Our GMRT survey provides 
systems in a narrow redshift range in which this measure- 
ment can be done. Thus high resolution optical spectroscopy 
of the corresponding QSOs are suitable to perform this test 
at 2; ~ 1.3. 
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